Identifying Nearby UHECR Accelerators using UHE (and VHE) Photons 
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Ultra-high energy photons (UHE, E > 10 19 eV) are inevitably produced during the propagation 
of ~10 20 eV protons in extragalactic space. Their short interaction lengths (<20 Mpc) at these 
energies, combined with the impressive sensitivity of the Pierre Auger Observatory detector to these 
particles, makes them an ideal probe of nearby ultra-high-energy cosmic ray (UHECR) sources. We 
here discuss the particular case of photons from a single nearby (within 30 Mpc) source in light of 
the possibility that such an object might be responsible for several of the UHECR events published 
by the Auger collaboration. We demonstrate that the photon signal accompanying a cluster of a 
few > 6 x 10 19 eV UHECRs from such a source should be detectable by Auger in the near future. 
The detection of these photons would also be a signature of a light composition of the UHECRs 
from the nearby source. 



The origin of ultra-high-energy cosmic rays (UHECR, 
E > 10 19 eV) still remains a mystery. However in recent 
ears clear progress has been made, with both the HiRes 
1] and the Pierre Auger Observatory (PAO) [2| provid- 
ing strong evidence of a suppression in the cosmic ray 
spectrum above ~ 5 x 10 19 eV. Moreover, the PAO has 
presented the first evidence that the arrival distribution 
of cosmic rays with energies larger than ~6x 10 19 eV is 
anisotropic [3| . If the flux suppression is interpreted as a 
consequence of severe proton energy losses during their 
propagation in extragalactic space, ie. the GZK feature 
0, [|, the sources of > 6 x 10 19 eV UHECRs would be 
required to reside in the nearby Universe, within several 
hundred Mpc of Earth. This anisotropy could therefore 
be a natural consequence of the fact that the local uni- 
verse is inhomogeneous. 

Presently indirect measurements of the UHECR com- 
position provide an unclear picture. Measurements of the 
AT m ax value of air showers detected by the PAO Q seem 
to indicate an increasingly heavy composition of UHECR 
whereas measurements by HiRes indicate a predomi- 
nantly light composition. Such measurements are of cru- 
cial importance in establishing whether the sources accel- 
erate predominantly protons or heavier nuclei and hence 
in determining the expected deflection of UHECR in both 
Galactic and extragalactic magnetic fields. Though the 
anisotropy signal observed by the PAO is perhaps more 
simply explained if the composition is light, interpreta- 
tions based on a few nearby sources of heavy nuclei are 
also possible 0. In either case propagation losses limit 



the UHECR horizon distance to below ~T00 Mpc at en- 
ergies >10 20 eV. In this letter we assume that UHECRs 
with energies > 6 x 10 19 eV are protons. 

For UHECR protons, the main channel of energy losses 
at energies >5x 10 19 eV is photopion production, with 
~0.5 mb cross-sections and inelasticities of the order 
m^/nip at threshold. Following their production and es- 
cape from the source region, UHECR protons undergo 
these photo-pion interactions with background radiation 
fields (predominantly the CMB) producing both neutral 
and charged pions and hence a secondary UHE photon 
and neutrino flux around the source region. With the 
propagation of these UHE-photons in turn being limited 
through 77 interactions with background photons, a re- 
gion of space should exist around the source in which 
this UHE-photon flux is maximal. With the 77 cross- 
section peaking at a value of roughly 200 mb, and the 
dominant photon background target to these photons, 
the cosmic radio background (CRB), having a density of 
~1 cm -3 , the interaction lengths for this process, above 
threshold, are also ~ a few Mpc. Thus, the distance for 
which the generated UHE-photon flux is maximum would 
lie between a few Mpc and a few tens of Mpc from the 
source. We note here that though the CMB presents a 
much larger photon background target in terms of num- 
ber density, the interaction of UHE-photons (and UHE- 
electrons) with it is suppressed due to the Klein Nishina 
suppression in the cross section for center-of-mass ener- 
gies ^> m e c 2 . 

For sources more than a few tens of Mpc away the 
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energy in the produced UHE photon population may 
go into the development of a full electromagnetic cas- 
cade, with the energy flux arriving mainly at GeV-TeV 
energies, such sources either contributing to the diffuse 
gamma-ray background at these energies, or being in- 
dividually resolvable in some scenarios @, [nj, depend- 
ing on the strength of the inter/extragalactic magnetic 
fields. However, for a more nearby source, the electro- 
magnetic cascade's development has insufficient time to 
develop, with only the appreciable production of a first 
generation of UHE photons, leading to their contribu- 
tion to the diffuse UHE photon background. Thus, the 
search for UHE-photons in the arriving UHECR flux can 
provide information about the local UHECR sources on 
these scales [H[. The level of this background at 10 19 eV 
is already strongly bounded by PAO observations to be 
less than ~2% of the proton flux [12J. This value ex- 
ceeds by a factor between a few and a few tens that ex- 
pected from theoretical calculations [3, 14 1. However, 
as demonstrated by 14J, an enhancement of the UHE 



ation rate of first generation of neurtral pions and the 
pair creation rate of UHE-photons, the distribution of 
UHE-photons with distance?, assuming rectilinear prop- 
agation, is, 



photon flux can be expected in the presence of a local 
source. 

A nearby source may easily result in the production 
of multiplets of UHECR events coming from within a 
few degrees on the sky, provided the composition at the 
source is dominated by protons at the highest energies 
to avoid sizable deflections in the intervening intergalac- 
tic magnetic fields (IGMFs). Furthermore, it has been 
suggested [l5[ that several of the 27 PAO events with 
energy > 6 x 10 19 eV used in the anisotropy study origi- 
nate from the nearest active galaxy - Centaurus A. At a 
distance of only ~3.8 Mpc [16| and with an exposure to 
the PAO of 1200 km 2 sr, if 2 of the PAO events originate 
from Cen A, the required luminosity to account for these 
events, assuming a source spectrum with spectral index 
2 between 10 19 eV and the cutoff, is ~ 0.4 x 10 40 erg s^ 1 . 
This is a small fraction of its total bolometric luminos- 



ity 12| and comparable to the inferred TeV luminosity 
as recently reported by the H.E.S.S. collaboration [l8|. 
In this way, Cen A may be considered as a promising 
UHECR source candidate. Indeed the possible subse- 
quent (cascade) TeV photon and UHE neutrino fluxes 
from Cen A have recently been calculated for this sce- 
nario 1!) . However, we note that it is still not clear if 
any region, within an object such as Cen A (jets, inner 
and outer lobes, black hole) meets the basic criteria for 
UHECR acceleration (see e.g. [13, [SH E3j ) . 

Following the analytic description of the photon flux 
development in for a very nearby (< 10 Mpc) ob- 
ject such as Cen A, a simple expression can be used 
to estimate the UHE photon flux coming from it. For 
this source most of the arriving protons have yet to un- 
dergo pion production interactions and the photon flux is 
dominated by the first generation of photons generated 
through 7T° decay. Through a consideration of the two 
coupled differential equations describing both the gener- 
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for which l pl is the interaction length for an UHE-proton 
of energy E p , ? 77 the interaction length for an UHE- 
photon of energy E 7 , the two energies being related 
through E-y « 0.5K P E P , where K p is the inelasticity of 
the j»7 collision [23j ]. In this expression 7V 7 (_B 7 ,Z) is the 
number of photons of energy £L at distance I from the 
source, similarly N p (E p ,0) is the number of protons of 
energy E p injected at the source. 

In order to more accurately describe proton propa- 
gation as well as explore the range of possibilities for 
the propagation of electromagnetic energy through in- 
tcrgalactic space we here employ a Monte-Carlo (MC) 
method for UHECR proton and photon/electron propa- 
gation as described in [3] . In this MC UHE photon pair- 
production interactions with background photons as well 
as UHE electron/positron production through charged 
pion decay are considered, both giving rise to a flux 
of secondary UHE-electrons. The environment through 
which these electrons propagate determines the ratio of 
the energy that these electrons input into a new gener- 
ation of UHE photons via inverse Compton (IC) scat- 
tering on CRB photons to that lost in much lower en- 
ergy synchrotron emission through interactions with the 
IGMFs. Both the strength of the IGMFs in the relevant 
volume and the density of the CRB are very uncertain. 
For our calculations, intergalactic magnetic fields in the 
range 10 -9 — 10~ 125 G are considered, (B-fields weaker 
than 10~ 125 G play no significant role for UHE-electron 
cooling on 30 Mpc distance scales). However, we note 
that the filling factor of IGMFs in galaxy cluster regions 
carry a large degree of uncertainty. We here assume that 
the sub /xG central field [24| drops off sufficiently quickly 
with distance that UHE-electron propagation occurs pre- 
dominantly outside such regions. We consider radio back- 
grounds ranging from the radio component of the CMB 
alone up to the measured, and possibly foreground con- 
taminated, value of [25[. For the range of radio back- 
ground values we consider, the uncertainty in the inter- 
action length of 10 19 eV photons is w25%, increasing by 
a factor of ^2 for photons with energy 3 x 10 19 eV. In the 
results that follow we account for the uncertainty in both 
the CRB and IGMFs, by showing the complete range of 
photon fluxes obtained from the ranges of both of these 
quantities. 

Fig. Q] shows the expected integrated rate (per year) 
of UHECR protons and photons from a source 3.8 Mpc 
away and 30 Mpc away arriving to the (complete) PAO. 
A power law injection spectrum at source with index 
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FIG. 1: The expected UHE photon arrival rate from a single 
source 3.8 Mpc (upper plot) and 30 Mpc (lower plot) away. 
For these plots, the arriving UHECR flux has been normalised 
such that 10% of the PAO events 0] above 6 x 10 19 eV came 
from the source. An injection index a — 2 and a cut-off energy 
E 

max — 10 20 5 eV were assumed. The uncertainty in the UHE 
photon flux is represented by the width of the band drawn for 
the photon flux. For the very nearby source, 3.8 Mpc away, 
the analytic result for the photon flux is shown for compari- 
son. These results bound the source distance range for which 
a local UHECR source leads to a strongly enhanced UHECR 
photon fraction, as demonstrated in Fig. 9 of [3] 



a =2 and an exponential cut-off at 10 20 5 eV was as- 
sumed. The flux of UHECRs from the source has been 
normalised such that it accounts for 10% (a doublet) of 
events above 6 x 10 19 eV recorded by the PAO in the 
data set Q (corresponding to just over 2 events in the 
data which had an exposure of 7000 km 2 steradian years). 
Note that the luminosity required by an UHECR source 
normalised in this way varies with declination, due to the 
non-uniform sky coverage of the PAO. In the first panel 
of Fig. [1] we show that this method gives a very similar 
result to the analytical description given above. 

Fig. [2] shows the impact on the expected arriving pho- 
ton rate from a source 3.8 Mpc away for different as- 
sumed source spectra. In the upper panel different injec- 



FIG. 2: The expected protons and photon arrival rates as 
for Fig. [T] from a source 3.8 Mpc away for which both the 
injection cut-off energy (£ max ) and injection index (a) have 
been varied. 



tion spectral indices of a — 2.4 and a — 2.7 are shown 
for comparison. In the case of a steep injection spec- 
trum produced by the source, with a — 2.7, a large con- 
tribution to the arriving UHECR flux at lower energies 
is expected. Such a large contribution from a nearby 
source (« 10%) over a large energy range would likely 
also lead to a detectable degree of anisotropy in the UHE- 
CRs with energies < 6 x 10 19 eV. Furthermore such a 
scenario would place uncomfortable requirements on the 
required source luminosity unless the spectrum becomes 
harder (flatter) at lower energies. As expected, the effect 
of a different injection spectrum on the photon rate is 
weak since the rate of UHECR protons have here been 
normalised to the integral flux at energies just above the 
photo-pion threshold and hence UHE photon production. 
The effect of reducing the maximum acceleration energy 
in the source, however, results in a significant reduction 
in the photon flux, as can be seen in the lower plot of 

Fig. m 

Accompanying this flux of UHE photons arriving to 
the Galaxy, a secondary flux of UHE electrons is also 
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generated. This is similar in magnitude to the UHE pho- 
ton flux, though it cuts off at energies somewhat below 
the equivalent cut-off in the UHE photon flux. These 
UHE electrons have similar (within a factor of 2) inter- 
action lengths compared to those of the UHE photons 
of the same energy as they propagate through the inter- 
galactic medium between Cen A and our Galaxy. Being 
charged particles these electrons will undergo deflections 
in the IGMFs. For - 1CT 10 G IGMFs, as suggested by 
[jjj ]. ~ 0.5° size deflections of 10 19 eV electrons can be 
expected for propagation distances of ^Mpc (assuming 
the field is coherent on this scale). 

Once the UHE electrons have entered the magnetised 
halo of our Galaxy, they can be expected to lose en- 
ergy through synchrotron emission in a time ti oss ~ 
1.3 B~ 2 Eig years, where is the magnetic field in 
/LtG and E\g is the electron energy in units of 10 19 eV. 
All of the electron's energy is lost in a small fraction, 
within 0.002°, of its gyration period, rapidly converting 
its energy into ~ 1.5 B^E\ 9 TeV photons. It is there- 
fore possible that the detection of UHE photons from a 
nearby source may be accompanied by a corresponding 
TeV gamma-ray flux originating from the same direction. 
The angular extension of this emission most likely being 
dominated by the angular spread of the arriving UHE 
electrons, expected to be ~ 0.5° for 10~ 10 G IGMFs. 
This 7-ray signal may be detectable by next generation 
Cherenkov telescope arrays such as CTA [27J ■ 

From figure [U the expected rate of > 10 19 eV photons 
from a source 3.8 Mpc away, such as Cen A, is seen to 
be 0.2 — 0.3 events/year in Auger assuming it is respon- 
sible for 10% of the > 6 x 10 19 eV flux. Should a large 
fraction of events in the current PAO dataset originate 
from Cen A, this implied rate will of course increase ac- 
cordingly. However, the discrimination of photons from 
protons does not have 100% efficiency. If the data from 
fluorescence detectors is used the efficiency is reduced 
by an order of magnitude due to their ~10% duty cy- 
cle. Surface detectors have a 100% duty-cycle but the 
overall photon efficiency (for a background free sample) 
in current analyses is close to 25% [28[ . No UHE photon 
candidate from anywhere on the sky has so far been iden- 
tified. In the search for these photons from an individual 
source their selection criteria could be significantly re- 
laxed - the background from misidentified protons should 
be less than ~1% in an 0.1 steradian region around the 
source - rather than in the 7r steradian field-of-view of 
Auger. A photon signal from single nearby sources is 
therefore potentially detectable using the PAO with a 
few years of data. 

Given the small statistics involved it is unlikely that 
the non- detection of UHE photons from an object like 
Cen A would place strong constraints on UHECR origin 
and composition in the near future. However, the detec- 
tion of one or more > 10 19 eV photons from a nearby ob- 
ject, particularly in combination with ~ 10 20 eV UHECR 



correlated to it, would provide extremely strong evidence 
for that object being an UHECR proton source. Along 
with this, the angular separation of such spatially cor- 
related UHE photon and UHECR events could provide 
constraints on the intervening Galactic and IGMFs. Fur- 
thermore, should future observations detect the associ- 
ated TeV emission from an UHE photon and UHECR 
source, the angular spread of this TeV emission may also 
allow the determination of the strength of intervening 
IGMFs. 

Given the importance of UHE photons as a diagnos- 
tic tool in UHECR studies we believe that effective dis- 
crimination between photons and nuclei should be a ma- 
jor design criterion for future UHECR detectors such as 
JEM-EUSO m and Auger North 0. 
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